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Abstract 

In [Schuhmacher, Electron. J. Probab. 10 (2005), 165-201] estimates of the Barbour- 
Brown distance d2 between the distribution of a thinned point process and the distribution 
of a Poisson process were derived by combining discretization with a result based on Stein's 
method. In the present article we concentrate on point processes that have a density with 
respect to a Poisson process. For such processes we can apply a corresponding result directly 
without the detour of discretization and thus obtain better and more natural bounds not 
only in d2 but also in the stronger total variation metric. We give applications for thinning by 
covering with an independent Boolean model and "Matern type I" -thinning of fairly general 
point processes. These applications give new insight into the respective models, and either 
generalize or improve earlier results. 



1 Introduction 

We consider thinnings of simple point processes on a general compact metric space X, where 
simple means that the probability of having multiple points at the same location is zero. The 
thinning of such a process ^ according to a [0, l]-valued measurable random field vr on Af is the 
point process ^tt, unique with regard to its distribution, that can be obtained in the following 
way: for any realizations ^{u>) (a point measure on M.^) and 7r{uj, •) (a function M.^ — > [0, 1]), look 
at each point s of ^(w) in turn, and retain it with probability 7r{u!, s), or delete it with probability 
1 — 7r{uj, s), independently of any retention/deletion decisions of other points. Regard the points 
left over by this procedure as a realization of the thinned point process ^j^. For a formal definition 
see Section [2j We usually refer to as the original process, and to tt as the retention field. 

The following is a well-established fact: if we thin more and more, in the sense that we con- 
sider a sequence of retention fields ('/r„)„gN with sup^g_;^' vr„(x) — > as n — > oo, and compensate 
for the thinning by choosing point processes ^„ whose intensity increases as n goes to infinity in 
a way that is compatible with the retention fields, then we obtain convergence towards a Cox 
process. The theorem below was shown in [12] for constant deterministic iTn, and generalized in 
f5] to general vr^. In order to specify what choice of the sequence (^n) is compatible with (vr^), 
we introduce the random measure A„ that is given by An{A) := J^7r„(x) ^n{dx) for every Borel 
set A in X. Convergence of random measures, and in particular of point processes, is defined via 
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the convergence of expectations of bounded continuous functions, where continuity is in terms 
of the vague topology (for details see [13j, Section 4.1). 

Theorem l.A (Kallenberg, Brown). For the sequences {iT„)nef>i '^''^d (^n)neN introduced above, 
we obtain convergence in distribution of the thinned sequence ((Cn)7r„)„gi^ towards a point process 

T] if and only if there is a random measure A on X such that A„ A as n ^ oo. In this case 
rj ~ Cox(A), i.e. rj is a Cox process with directing measure A. 

In [18] the above setting was considered for the situation that the are point processes 
on [0, 1]^ that are obtained from a single point process ^ by gradual contraction of using 
the functions k„ given by Hn{x) ■= {l/n)x for every x € (for notational convenience in 
the proofs, the order of contracting and thinning was interchanged). Under the additional 
assumption that ^ and 7r„ satisfy mixing conditions, which makes it plausible for the limiting 
process in Theorem II. Al to be Poisson (see the remark after Theorem l.A of [E]), several 
upper bounds for the Wasserstein-type distance d2 between the distribution of the thinned 
and contracted process and a suitable Poisson process distribution were obtained under various 
conditions. These results were derived by discretizing the thinned process and the limiting 
Poisson process, and applying then a discrete version (essentially formulated as Theorem 10. F 
in [1]) of the "local Barbour-Brown theorem", Theorem 3.6 in which is based on Stein's 
method (see [20] and [3]). 

Although the bounds were of good quality and have proved their usefulness in several appli- 
cations, they had some shortcomings, which were mainly related to the fact that they were still 
expressed in terms of discretization cuboids. This made the results rather unpleasant to apply 
in many situations where truly non-discrete point processes were considered. 

The present article deals with these issues by offering a proof that makes direct use of 
Theorem 3.6 in p] without the intermediate step of discretization. There are several advantages 
of this approach: the proofs become simpler and more elegant; the upper bounds are much 
more natural, more easily applied, hold for the most part under more general conditions and are 
qualitatively slightly better; what is more, we can obtain a bound not only for the d2-distance, 
but also for the stronger total variation distance. In order to apply Theorem 3.6, we restrict 
ourselves to point processes ^ that have a density with respect to the distribution of a simple 
Poisson process, which is a natural and common choice for a reference measure and leaves us 
with a very rich class of processes. The simplicity of the Poisson process implies the simplicity 
of ^. It should be noted, however, that any non-simple point process can be turned into a simple 
one by lifting it to a bigger space X x K. This space has to be compact in our setting, so that 
e.g. K = {1,2, . . . , A;} (if there is a maximal number k of points per location), K = [0, 1], or 
K = [0, e] for e small are typical choices. 

We start out in Section [2] by giving an overview of the technical background and some 
notation needed to formulate and prove the main results. These results are then presented 
in Section 3. We provide upper bounds for the dxv- and the d2-distances between ^{^n) 
and a suitable Poisson process distribution, first in a general setting, and then for a number 
of important special cases. The last of these special cases (see Corollary 13. Ep is suitable for 
comparison with the upper bounds in [18]. Finally, in Section U two "extreme" applications of 
the main results are studied, in which the retention field takes only the values and q E [0, 1] and 
is either independent or completely determined by ^. In the first application a point is always 
deleted if covered by a certain independent Boolean model, and is retained with probability q 
otherwise. In the second application, a point is always deleted if there is any other point present 
within a fixed distance r (following the construction of the Matern type I hard core process), 
and again is retained with probability q otherwise. 
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2 Preliminaries 



We first introduce some basic notation and conventions, before giving an overview of some of the 
theoretical background and presenting the more technical definitions in the various subsections. 

Let {X, do) always be a compact metric space with do < 1 that admits a finite diffuse measure 
a 7^ 0, where diffuse means that a{{x}) = for every x G A". Denote by B the Borel <T-algebra 
on X, and by Ba the trace cr-algebra B\a = {B H A; B £ B} for any set A C X. Furthermore, 
write 971 for the space of finite measures on X, and equip it with the vague topology (see [13], 
Section 15.7) and the corresponding Borel fj-algebra M (see [131 . Lemma 4.1 and Section 1.1). 
Do the same for the subspace C 9Jt of finite point measures, and denote the corresponding cr- 
algebra by A/". Write furthermore := {g G 9T; q{{x}) < 1 for all x € X} for the A/"- measurable 
set of simple point measures. A random measure on A:' is a random element of SOT, and a point 
process on X a random element of A point process ^ is called simple if G = 1. By 
Po(A) we denote the distribution of the Poisson process on X with intensity measure A if A G Tl, 
and the Poisson distribution with parameter A if A is a positive real number. 

We think of measures (random or otherwise) always as being defined on all of X. Thus, 
for any measure fj, on X and any A B, we denote by fi\A the measure on X that is given 
by ij.\a{B) := n{A n B) for all B e B. Let Tl{A) := {/uU; /i G 9?t}, ^{A) := {^U; Q G OT}, 
M{A) := M\m{A) = {C n m{A); C G M}, and M{A) := M\^^a) = {C n m{A)- C G M}. 
Furthermore, set ^*{A) := ^{A) n9T*. Sometimes absolute value bars are used to denote the 
total mass of a measure, i.e. |^| := n{X) for any G dJt. 

For 0" G 91*, we do not notationally distinguish between the point measure and its support. 
Like that we can avoid having to enumerate the points of the measure, which sometimes saves 
us from tedious notation. Thus, the notation J f da = X^seo- /('^) used instead of writing 

J f da = Yll=i /(^i)) where a = J2^=i ^s^- The same concept is extended to the realizations 
of a simple point process, as long as it is not important what happens on the null set of point 
patterns with multiple points. Hence we may also write E(/ f d^) = E{T,sei; /(«)) if ? is simple. 
In order to facilitate the reading of certain formulae, we furthermore make the convention of 
using the letters x,x,y for general elements of the state space X, while s,s,t are reserved for 
points of a point pattern in X. Finally, we sometimes omit the addition "almost surely" or 
"(for) almost every ..." for equations and inequalities between functions on measure spaces if 
it is evident and of no importance that the relation does not hold pointwise. 

2.1 Densities with respect to the standard Poisson process distribution Pi 

Let a 7^ be a fixed diffuse measure in dJt, which we will regard as our reference measure 
on X. Typically, if (a superset of) X has a suitable group structure, a is chosen to be (the 
restriction of) the Haar measure. If X C M^, we usually choose a = Leh^\x of course. We 
write Pi := Po(q) for the distribution of what we call the standard Poisson process on X, and 
Pi^A ■= Po(aU) for the distribution of the Poisson process on X with expectation measure a\A- 
It is convenient to admit also a{A) = 0, in which case -Pi, a = ^O) where denotes the zero 
measure on X. 

A popular way of specifying a point process distribution is by giving its Radon-Nikodym 
density with respect to the distribution of the standard Poisson process (if the density exists; 
see p[5. Sections 6.1 and 6.2 for a number of examples). The following simple lemma will be 
useful. 

Lemma 2. A. For any a > 0, a density /a of Pa^A '■= Po(a-a|A) with respect to Pi, a is given by 



3 



for every g G ^{A) . 

Proof. See [H], Proposition 3.8(ii), for the case X C M^, a = Leb^lx, and a more general 
Poisson process. The same proof holds for general X and a. □ 

To avoid certain technical problems, we require our densities to be hereditary whenever we 
are dealing with conditional densities. 

Definition. A function / : 91 — > M+ is called hereditary if f{Q) = implies f{cr) = whenever 
^, CT G with g C a. 

The point processes on X that have a hereditary density with respect to Pi form the class of 
Gibbs point processes. These include pairwise and higher order interaction processes (see [H], 
Section 6.2). The general form of a Gibbs process density is given in Definition 4.2 of [1]. 



2.2 Thinnings 

In what follows, let ^ always be a point process on X that has density / with respect to Pi, and 
let TT := {tt{-,x);x G <Y) be a [0, 1] -valued random field on X that is measurable in the sense 
that the mapping Q x X ^ [0,1], {to, x) i-^ 7r(a;, x) is (7{tt) (g) ;B-measurable. In the main part we 
strengthen the technical conditions on ^ and vr somewhat in order to avoid some tedious detours 
in the proofs. 

We use the definition from [18j for the vr-thinning of ^. 

Definition (Thinning). First, assume that ^ = o" = and tt = p are non-random, where 

V G Si G X, and p is a measurable function X ^ [0; !]• Then, a vr-thinning of ^ is defined as 
= Yli=i-^i^Si, where the Xi are independent indicator random variables with expectations 
p{si), respectively. Under these circumstances, has a distribution P(cr,p) that does not depend 
on the chosen enumeration of a. We obtain the general vr-thinning from this by randomization, 
that is by the condition P[^^ G • | '^l = ^(5,7r) (it is straightforward to see that P((^n) is a 
(t(^, vr)-measurable family in the sense that P(^^^^^{D) is (t(^, vr)-measurable for every D G 9T). 
Note that the distribution of is uniquely determined by this procedure. 

The following lemma gives the first two factorial moment measures of ^t^. For g = Y2J=i £ 
DT, write fjt^l := Yll j=i i^j ^is^,sj) ^or the factorial product measure on X x X. Remember that 
the expectation measure fii of ^ is defined by fJ-iiA) := E(^(^)) for every A £ B, and that the 
second factorial moment measure /Upj of ^ is defined to be the expectation measure of^I^l. Let 
A be the random measure on X that is given by A(^) := J^vr(x) ^{dx) for A £ B (cf. A„ in 
Section [T|). 

Lemma 2.B. We obtain for the expectation measure /i^^' and the second factorial moment 
measure ^^j' o/^,r 

(i) fi^^\A) =^(^j ^(a;) ^{dx)^ = E(A(^)) for every A G B; 

(ii) /ujJ(S) = E^^y vr(x)vr(x) i^^\d{x,x))^ =: E(aP1(P)) for every B G S^. 

Proof. Write ^ = ^^=1 ^Sn where V and Si are cr(i^)-measurable random elements with values 
in and X, respectively. Such a representation exists by Lemma 2.3 in [13j. 
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(i) For every A B we have 



(ii) For every B G we have 

=E(d'](i?)) =E(E(ei2](i?) I C,^)) 

□ 

2.3 Metrics used on the space of point process distributions 

We use two metrics on the space of probabihty measures on The one that is more 

widely known is the total variation metric, which can be defined on any space of probability 
measures. For P, Q G *p(9T) it is given by 



or, equivalently, by 



dTviP,Q) = sup|P(C)-Q(C)| 

CGAf 

dTy(P,g) = min P[6/6]- (2.1) 

See [1], Appendix A.l, for this and other general results about the total variation metric. 

The second metric we use is a Wasserstein type of metric introduced by Barbour and Brown 
in [2], and denoted by d2- It is often a more natural metric to use on *p(9T) than d-ry, because 
it takes the metric do on X into account and metrizes convergence in distribution of point 
processes. The total variation metric on the other hand is strictly stronger, and at times too 
strong to be useful. 

Define the di-distance between two point measures = Yli=i i and Q2 = X]l=l <^S2 i ^ 

as 

'l if Iftl 7^ \Q2\, 

di{Qi,g2) ■■= < min^es„ [^Ei'=Wo(si,i,S2,r(i))] if l^*!! = l^'sl = v > 0, (2.2) 
,0 if |£>i| = \g2\ = 0, 

where denotes the set of permutations of {1,2, ... ,v}. It can be seen that (OT, di) is a 
complete, separable metric space and that di is bounded by 1. Now let ^2 := {/ : ^ — > 
^'1 \f{Qi) ~ fiQ2)\ < di{gi, Q2) for all ^1, Q2 S 01}, and define the d2-distance between two mea- 
sures P,Q e *P(0T) as 

d2{P,Q) := sup [ fdP- [ fdQ 
fdTi J J 

or, equivalently, as 

d2(P,g) = min Edi(ei,6)- (2.3) 

See |18j, [19j . and [22] for this and many other results about the Barbour-Brown metric ^2- By 
()2.ip and ()2.3p we obtain that d2 < dry. 
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2.4 Distance estimates for Poisson process approximation of point processes 
with a spatial dependence structure 

In this subsection, a theorem is presented that provides upper bounds for the total variation and 
d2 distances between the distribution of a point process with a spatial dependence structure and 
a suitable Poisson process distribution. This result is very similar to Theorems 2.4 and 3.6 in [2], 
but deviates in several minor aspects, two of which are more pronounced: first, we use densities 
with respect to a Poisson process distribution instead of Janossy densities, which simplifies 
certain notations considerably; secondly, we take a somewhat different approach for controlling 
the long range dependence within ^ (see the term for (j){x) in Equation ()2.1ip ). which avoids the 
imposition of an unwelcome technical condition (cf. Remark I A. Cp . 

Let ^ be a point process on X whose distribution has a density / with respect to Pi and 
whose expectation measure n = ni is finite. Then ^ has a density u with respect to a that is 
given by 

u{x)= [ f{Q + 6^) Piidg) (2.4) 

for a-almost every x (z X, which is obtained as a special case of Equation ()2.8p below by choosing 
AT, := {x}. 

For A ^ B, we set 

fA{Q) := [ f{Q+Q)PiMdQ) (2-5) 

for every g G ^(A), which gives a density : ^ M+ of the distribution of with respect 
to Pi^A (we extend /^i on \ ^{A) by setting it to zero). This can be seen by the fact that 
integrating /a over an arbitrary set C G ^(A) yields 

/ fA{Q) PiAdQ) = [ I G C]f{g + g) PiMdg) PiAido) 

Jc J'niA) j?n{A=) 

= [ l[a\A G C]f{a) Piida) 

where we used that (^U;??U=) ^ ® Pi,A'' for r/ ~ Pi. Note that the argument remains 
correct if either a{A) or a{A'^) is zero. 

We introduce a neighborhood structure {Nx)xex on X, by which we mean any collection of 
sets that satisfy x E for every x ^ X (note that we do not assume Nx to be a neighborhood 
of X in the topological sense). We require this neighborhood structure to be measurable in the 
sense that 

N{X) := {(x, y) ^ X^;y e iVj G B\ (2.6) 

This measurability condition is slightly stronger than the ones required in [2] (see the discussion 
before Remark 2.1 in [7J for details), but quite a bit more convenient. The play the role of 
regions of strong dependence: it is advantageous in view of Theorem 12. CI below to choose 
not too large, but in such a way that the point process ^ around the location x depends only 
weakly on ^|7v=- Write Nx for Nx \ {x}. 

We use the following crucial formula about point process densities, which is proved as 
Proposition lA.AI in the appendix. For any non-negative or bounded measurable function h : 
X X M, we have 

El/ h{x,(,\Ng) C{dx)] = / h{x,g)fNc(j{x}{Q + Sx) Pi,Niidg) a{dx). (2.7) 
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As an important consequence we obtain by choosing h(x, g) := IaIx) that 

u{x)=l fmu{x}iQ + Sx) Pi,Ns{dg) (2.8) 

for a-almost every x € X. 

In many of the more concrete models, the density / of ^ is hereditary and therefore ^ is a 
Gibbs process, in which case the above expressions can be simphfied by introducing conditional 
densities. Let K := Uxexil^} ^ ^i^x)) C ^ x 91 and define a mapping 51 : ^ x — > by 

9{x; q) := ) \ 2.9 

for {x, g) € ^ and g{x; g) := otherwise, where the fraction in ()2.9p is taken to be zero if the 
denominator (and hence by hereditarity also the enumerator) is zero. For (x, g) (z R the term 
g{x; g) can be interpreted as the conditional density of a point at x given the configuration of 
^ outside of is g. Equation ()2.7p can then be replaced by the following result, which is a 
generalization of the Nguyen-Zessin Formula (see [TB], Equation (3.3)): for any non-negative or 
bounded measurable function /i : <Y x 91 — M, we have 

E(^j^h{x,^\N^J ^dx)^ = j^E{h{x,C\NM^-^^\m)) a{dx). (2.10) 

This formula was already stated as Equation (2.7) in [2j for functions h that are constant in x and 
as Equation (2.10) in [7] for general functions, both times however under too wide conditions. 
See Corollary lA.BI for the proof and Remark lA.CI for an example that shows that an additional 
assumption, such as hereditarity, is required. As an analog to (|2.8|) . we obtain 

u{x)=E{gix;C\Ns)) 

for a-almost every x ^ X. 

We are now in a position to derive the required distance bounds. 

Theorem 2.C (based on Barbour and Brown [2], Theorems 2.4 and 3.6). Suppose that ^ is a 
point process which has density f with respect to Pi and finite expectation measure /i. Let fur- 
thermore {Nx)xex be a neighborhood structure that is measurable in the sense of Condition l\2.Q\] . 
Then 

(i) dry(^(0,Po(/x)) < j^ii{N,) ^x{dx)+¥.l^j^i{Nx)i{dx)^ + j J{x) a{dx); 

+ Mi{n) I 0(x) a{dx); 



li{N^) ^^{dx) + ^ / i{N^)i{dx] 
X \Jx 



(11) d2(^(0,Po(/x)) < M2(/i) 

where 

Mi(/x) =min(^l,^j^^, M2(/i) = min (^1,^(^1 + 2 log+ 



X 
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and 



4^{x) = I |/iV|u{x'}(^ + ^x) - fNg{Q)u{x)\ Pi,N^{dg) 



2 sup 

CeM{Ni) 



c 



[fNgu{x}{Q + ^x) - fN^{Q)u{x)] Pi,Ns{dg) 



(2.11) 



If f is hereditary, (p can be rewritten as 

4>{x)=E\g{x;^\Ns)-u{x)\. (2.12) 
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Remark 2.D. We refer to the three summands in the upper bound of Theorem 2.C.(i) as basic 
term, strong dependence term, and weak dependence term, respectively. The basic term depends 
only on the first order properties of ^ and on a{Nx)- The strong dependence term controls what 
happens within the neighborhoods of strong dependence and is small if a{Nx) is not too big and 
if there is not too much positive short range dependence within Finally, the weak dependence 
term is small if there is only little long range dependence. 

Remark 2.E. Theorem 5.27 in [22] (which is based on several earlier results by various authors) 
gives an alternative upper bound for the (i2-distance above, which when carefully further esti- 
mated is in many situations superior to the one in [2], insofar as the logarithmic term in M2(/i) 
can often be disposed of. After applying the same modifications as in the proof of Theorem 12. CI 
below, it can be seen that 



d2(^(e),Po(;u)) 



< 



E 



+ E 



3.5 



+ 



2.5 



/"I am + 1 

3.5 2.5 



+ Mi(//) / (l){x) a{dx). 
'x 



Since working with this inequality requires a more specialized treatment of the thinnings in 
our main result, and since the benefit of removing the logarithmic term is negligible for most 
practical purposes, we do not use this bound in the present article. 



Proof of Theorem \2.CX Following the proof of Theorems 2.4 and 3.6 in [2] (applying Equa- 
tions (2.9) and (2.10) of [2j, but not Equation (2.11)), we obtain by using Stein's method 
that 



|E/(e)-E/(r?)| < As/i 
+ 



n{N^) ^i{dx) +^ / i{N^)i{dx] 

X \J X 



E / [/i(eU,^ + 5.) - h{i\N^)\ {i{dx) - fi{dx)) 
Jx 



(2.13) 



for every / G J^tv = {Ic; C G M} [or / G .7^2 in the case of statement (ii)], where r] ~ Po(/i) 
and h := are the solutions of the so-called Stein equation (see [2], Equation (2.2)), which 
have maximal first and second differences 



Aih := sup \h{g + Sx) — h{g)\ 
Qem,xex 



and 



A2/1 := sup \h{Q + dx + Sy) - h{g + 6x) - h{g + 6y) + h{g) 

g&%x,y<^X 



that are bounded by 1 [or Aih < Mi(fi) and < M2(/i) in the case of statement (ii); see 
Propositions 5.16 and 5.17]. 

All that is left to do is bounding the term in the second line of (I2.13p . which is done as 
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follows. Setting C+ := {g G 9T(iV^); fNiu{x}ig + ^x) > fNs{g)u{x)} G AA(A^^), we obtain 
IE [_[mNs + Sx) - mN^)] i^dx) - Kdx)) 

[h{g + Sx) - h{g)] {fN^u{x}{g + ^x) - fNg{g)u{x)) Pi,Ng{dg) a{dx] 



X 




X J'yi{N-) 

< AJi / / \fmu{x}{g + Sx) - fNiig)u{x) \ Pi,N^{dg) a{dx) 
Jx J'niNg) 

= 2Aih / {fN^u{x}{g + Sx) - fNi{g)u{x)) Pi,Ng[dg) a{dx) 
Jx Jet 



2 Aih I sup 

IX C<^{N^) 



c 



{fNsu{x}ig + Sx) - fNgig)u{x)) Pi,Ng{dg) 



a{dx) 



(2.14) 



where we use Equation ()2.7p for the second line and 

{fNsu{x}{g + ^x) - fN^{g)u{x)) Pi,N^{dg) = 



for the fourth line, which follows from Equation ([2.80 . The integrands with respect to a{dx) in 
the last three lines of (j2.14p are all equal to 



\fmvj{x}{'y\m +^^) - fm{'^\Niju{x)\ Pi{da) 



and hence S-measurable by the defintion of and the fact that Condition ()2.6p implies the 
measurability of the mapping [r^" x 91 — > x OT, (x, a) i— > (x,(T|Arc)] (see |73 , after Equation (2.4)). 

Pluging (I2.14|) into ()2.13p and taking the supremum over / completes the proof for general /. 
If / is hereditary, we have f n^\j{x}{q ~^ ^x) = g{x] g)fN^{g): so that the additional representation 
of claimed in (I2.12p follows from the third line of Inequality ()2.14p . □ 



3 The distance bounds 



We begin this section by presenting the general setting for our main result. Theorem 13. At partly 
compiling assumptions that were already mentioned, partly adding more specialized notation 
and conditions. Thereafter the main result and a number of corollaries are stated, and in the 
last subsection the corresponding proofs are given. 



3.1 Setting for Theorem [3TA] 

Let ^ be a point process on the compact metric space {X^do) which has density / with respect 
to Pi and finite expectation measure = ^i. Furthermore, let vr = {t[{-,x); x € Af) be a 
[0, l]-valued random field. We assume that vr when interpreted as a random function on X takes 
values in a space E C [0, 1]'^ which is what we call an evaluable path space (i.e. <^ : E y. X ^ 
[0,1], (p,x) I— > p{x) is measurable) and that there exists a regular conditional distribution of 
vr given the value of ^. Neither of these assumptions presents a serious obstacle; we refer to 
Appendix IA.3I for details. Let then A be the random measure given by A(^) := J^vr(3;) ^{dx) 
for ^ G i3, and set hP'^B) := /^vr(x)vr(x) Sp'^{d{x,x)) for any 5 G 0^ 

Choose a neighborhood structure {Nx)x(^x that is measurable in the sense of Condition (|2.6p . 
We assume for every x ^ X that vr(x) and t^\n^ are both strictly locally dependent on ^ in such 
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a way that the corresponding "regions of dependence" do not interfere with one another. More 
exactly, this means the fohowing: introduce an arbitrary metric do on X that generates the 
same topology as do, and write B(x,r) for the closed dg-ball at x G with radius r > and 
M{A,r) := {y £ X; do{y,x) < r for some x £ A} for the do-halo set of distance r > around 
A d X. Suppose then that we can fix a real number R > such that 

B(x, R) n ]B(iV^, i?) = (i.e. D M{x, 2R)) (3.1) 

and 

^5Ib(.,«) ^lB(a;,i?)= and t:\nc ^^i.^^.,^, i\nm,RY (3-2) 

for every x ^ X, where X JL^ Y denotes conditional independence of X and Y given Z. If Z is 
almost surely constant, this is just the (unconditional) independence of X and Y; in particular 
we require 7r(x) JL ^ if a{E{x,R)) = 0. Set ^int := A\^t{x) := B(x,i?) and ^ext := -^extix) ■ = 
B(A'^, R), where we usually suppress the location x when it is clear from the context. 

We introduce two functions to control the dependences in (^,vr). The function (3 : X ^ 
is given by 



j3{x) := / E(7r(x) I ^Ui„t = ftnt + ^x) / \f{QexU Q\nt + 5x)\ Pl,A,AdQeyit) PlAlnA'^Qh 

(3.3) 

where f{Qext-,Q\at + ^x) '■= /AcxtUAint(fet + ) — /Acxt(^'ext)/Aint(ffint + ^x), and hence 

controls the long range dependence within ^, as well as the short range dependence of vr on ^. 
If a{Aiat) = 0, the conditional expectation above is to be interpreted as E(7r(x)) for every 
^'int £ ^(^int)- The function 7 : ^ — > is taken to be a measurable function that satisfies 

/ esssup cov{TT{x),lc\i = Q + 5x) f{Q + 5x) Pi{dQ) <^{x), (3.4) 

and hence controls the average long range dependence within vr given ^. For the definition 
of the essential supremum of an arbitrary set of measurable functions (above, the functions 
i-^' |cov(7r(x), Ic = ^ + '5x)|] for C € (7(7r|iv^)), see [15], Proposition II. 4.1. 

3.2 Results 

We are now in the position to formulate our main theorem. 

Theorem 3. A. Suppose that the assumptions of Subsection \3.1\ hold and write /i*'^' = EA for 
the expectation measure of We then have 

(i) dTv{^{Cn),Po{fJ'^"^)) < {EAf{N{X)) +KA^^\n{X)) + [ (3{x)a{dx)+2 [ j{x) a{dx); 

Jx Jx 

(it) d2(^(eO,Po(^'"0) <^^2(/i<"')((lEA)'(iV(A'))+EA[- 



where 



+ Mi(/i'^))( / p{x)aidx) + 2 / j{x)a{dx) 
'x Jx 



1 Q^lI \ 

Mi(^'''^) = min( 1, " ) , and 



yJ^MX))'' 

M2(/.(-0 = min(l,^^^(l + 21og+(^EA(A')))), 
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If ^ and TT are independent, we obtain an interesting special case, where the upper bound can 
be expressed in terms of essentially the quantities appearing in Theorem I2.CI which are based 
solely on ^, and some rather straightforward quantities based on vr. 

Corollary 3.B (Independent case). Suppose that ^ is a point process on {X,do) which has 
density f with respect to Pi and finite expectation measure fi = fii. Let vr = (7r(-,x); x G 
be a [0, l]-valued random field that has an evaluable path space E C [0, 1]'^ and is independent 
of ^. Choose an arbitrary neighborhood structure {Nx)x€X (always x G Nx!) that is measurable 
in the sense of Condition (12. 6p , and take 7 : A" — > to be a measurable function that satisfies 

sup |cov(7r(a;), lc)| • u(x) < 7(x) (3-5) 

Ceo-(7r|ivc) 

for almost every x. Note that the expectation measure of^T^ is /U^'^'(-) = EA(-) = J_E7r(x) ^i[dx), 
and let Mi(^'-^^) and M2(//*'^') be defined as in Theorem \3.Al We then obtain 

(i) dTy(^(e.),Po(/i<"')) 

< / / E7r(x)E7r(x) /xi(dx) + / E(7r(x)7r(x)) ^upj x)) 
Jx Jn^ Jn{X) 

+ / 'Ett{x) (f>{x) a{dx) + 2 / ^{x) a{dx); 
Jx Jx 

(ii) d2(^(e.),Po(/i''^')) 

<M2(Ai*'")(/ / Ett{x)Ett{x) ixi{di) iii{dx) + E{tt{x)tt{x)) iim{d{x,x))\ 

\Jx Jn^ Jn(X) J 

+ Mi{fi^^^)(^j E7r{x)4>{x) a{dx) + 2j 7(3;) a{dx)^; 

where (p is given by Equation ()2.1ip and, if f is hereditary, by Equation ()2.12p . 

A further corollary is given for the case of a deterministic retention field, which means that 
the retention decisions are independent of each other given the point process ^. We formulate 
only a very special case, where the point process lives on R-^ and various spatial homogeneities 
are assumed, which leads to a particularly simple upper bound. The corollary also illustrates 
how we can deal with the issue of boundary effects in the neighborhood structure by extending 
the Nx beyond X in such a way that they are translated versions of each other and that the 
same holds true for their "complements" Mx \ Nx ■ 

In this article we always tacitly assume that is equipped with the Euclidean topology. 
Write \A\ := Leb^{A) for any Borel set A C R^, and X + X' := {x + x'; x £ X,x' e X'} and 
X — X := {x — X] X £ X} for X,X' C and x € M'^. For the definition of point processes 
on non-compact spaces and elementary concepts such as stationarity, we refer to the standard 
point process literature (e.g. [E] or [8]). 

Corollary 3.C (Constant case). Let X,y C M.^ be two compact sets, where X has positive 
volume and X + UsigA'C'^ — x) C y, and consider a metric dQ < 1 on X that generates the 
Euclidean topology (typically, dQ(x,y) = min(|x — y\, 1) for all x,y £ X). 

Suppose that ( is a stationary point process on whose restriction ^ := C\y is a Gibbs process 
with density f with respect to Pi = 'Po{Leh^\y) and finite expectation measure ji = ^i = 
miLeh^\y. Denote by /C the second reduced moment measure of C, (see around Equation (|4.7p 
for details). Let {Nx)x£X o,nd {Mx)x&x be two neighborhood structures whose sets Nx := x + N 
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and Mx := x + M (not necessarily C X now!) are translated copies of single hounded measurable 
sets N,M C which are chosen in such a way that N C M and X C Mx C 3^ for every 
X ^ X . Choosing our retention field ir = p £ [0, 1] to be deterministic and constant and noting 
that /x'^'Ia" = pmiLeh^lx , we then have 

(i) dTv{-^{ip\x),^o{pmiLeh^\x)) < p^ml\X\{\N\ + 1C{N)) + p\X\¥.\T 

(ii) d2(^(ep|A'),Po(pmiLeb^|;,)) < min(pmi|A'|, ^(1 + 2 log+ (^^^^)))pmi(|iV| + 1C{N)) 

+ min {^Jpmi\X\, 1.647) ^Jpmi\X\^^T - Er|, 

where T is an arbitrary random variable that has the same distribution as g' {x]^\m^\Nx) ~ 
fiM^\N^)u{x}{^\M^\N^ +4)//A/,\ArJ?|M,\7vJ (for One and therefore every x G X). 

Remark 3.D. While it seems very appealing to admit Mx = y = M^, this case actually re- 
quires a different and somewhat technically involved construction for the conditional density 
^|jgD\7Vi), because it cannot reasonably be assumed that a density of a point process dis- 
tribution with respect to the standard Poisson process distribution exists if the state space is 
(consider for example a hard core process: the hard core event that no two points are closer 
than some fixed distance r > is a Po(Leb^)-null set). As a matter of fact, the natural setting 
is that of a Gibbs process on the whole of M.^ defined via a stable potential on the set of finite 
point measures on M.^ , which essentially provides us with "compatible" conditional densities for 
the point process distribution on bounded Borel sets given the point process outside (see [T7] . 
from page 6.9 onwards, for a detailed construction). For a fixed bounded Borel set 3^ C we 
write /b(- | r) : 9T(i3) — > M+ for the conditional density of given ^\yc = r. It can then be 
seen that the crucial inequality 

^g{x;^\x\Nj - mi\ < E\g'{x;^\M^\Nj - 

(see Inequality (|3.14p ) and hence the proof of Corollary 13. CI can be reproduced under very 
general conditions if Mx = M.^ , where 



g{x;Q) = 

for every g € ^{X \ Nx), and 



I<n{yc) f{X\N^)u{x} ig + Sx\ t) ^{i\y<^) ^ {dr) 
f{y\N^)vj{x}{Q\y\N^ + I Q\y'=) 



g'{x; g) 



fy\NAe\y\N^ l^ly 



for every g G ^{W^\Nx) (as earlier we set such fractions to zero if the denominator is zero). By 
the construction of the Gibbs process on (using Equation (6.10) in [17]), the term g'{x;g) 
does not depend on the choice of y D Nx, except for ^3 in a P(.^|iRi5\jVi)~^-null set. It can be 
interpreted as the conditional density of a point at x given C|rd\7v = Q- 



In the next result, the situation of Theorem 13. Al and its corollaries is examined for the case 
where we compensate for the thinning by contracting the state space as it was done in |18j . 



Corollary 3.E (Thinning and contraction in R^). Suppose that X is a compact subset o/M^ 
and that a = Leh^\x- Let T > 1 and kt ■ — M^,x i— > ^x. Assume furthermore that the 
metric do on X generates the Euclidean topology and satisfies dQ(^KT{x), KT{y)) < do{x,y) for 
every x,y £ X. 
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Then Theorem \3.A\ and Corollaries \3.E\ and \3.Q remain true if the point processes on the left hand 
sides of the estimates are replaced by their respective image processes under the contraction kt . 
We thus obtain under the general prerequisites of Theorem \3.A\ 

(i) dTv{-^{i^K-T^).^o{ii^^^K^^)) < (EA)^(A^(;f)) + EAPI (Ar(;f )) + f p{x) dx + 2 [ ^{x) dx; 

Jx Jx 

(ii) d2(^(e.K^^),Po(^("'K^i)) <M2(^("')((EA)'(7V(A')) +EA[2l(iV(;f))) 

+ Mi{ii^^^)(^j dx + lj j{x)dxy 

Remark 3.F (Comparison with [18j). The setting of Corollary 13. El corresponds in large parts 
to the situation in [18], especially if we set X := k^^{J) for a fixed compact set J C and 
compare our statement (ii) to Theorem 3.B. 

It is more strict in essentially two respects. First, of course, we admit only point processes 
whose distributions are absolutely continuous with respect to a homogeneous Poisson process. 
Secondly, we require strict local dependence of vr on ^ (see Condition ()3.2p ). which in [18] was 
only done for Section 4 (in slightly different form), but which also gives the direct benefit of a 
conceptually simpler and more intuitive control of the long range dependences. 
On the other hand, the setting of Corollary 13. El gives us more freedom than we had in [18] in 
the sense that the objects live on a general compact subset of MP , that there are only minimal 
conditions on the moment measures (as opposed to Assumption 1 in [IB]), and that the choice 
of do and of the neighborhoods of strong dependence is much wider. 

Regarding the upper bounds obtained we have clearly improved. The terms in our statement (ii) 
above have their counterparts in the various terms in Theorem 3.B of [18] (with the integrals 
over f3{x) and 7(a;) being summerized as a single long range dependence term), but have become 
simpler and more natural, without any suprema or infima over discretization cuboids and with 
explicit and manageable constants. The bound as a whole is somewhat better (if we use a 
heuristic approach for comparing the long range dependence terms) and quite a bit more easily 
applied, which can be seen from comparing the application in Subsection 14.11 with the same 
application in Subsection 3.3 of [T8] . 

Remark 3.G. As pointed out in [T8], it may be desirable to approximate the distribution 
of the thinned point process by a Poisson process distribution that has a somewhat different 
expectation measure. Corresponding distance estimates can easily be obtained from upper 
bounds for distances between Poisson process distributions. We have cZtv (Po(A), Po(;u)) < 
||A — ^11 for A, G 571 by Remark 2.9 in [2], where ||-|| denotes the total variation norm for signed 
measures. For d2 an upper bound is given as Inequality (A. 3) in (which is the same as 
Inequality (2.8) in [6j). 

3.3 Proofs 

Proof of Theorem \3.Al By Lemma lA.DI in the appendix a density /^'^^ of ^-^ with respect to Pi 
exists, and the finiteness of /i implies the finiteness of /x^'^'. Hence we can apply Theorem 1 2. CI 

The integrals in the upper bound can be further evaluated as follows. For the first two 
integrals (basic term and strong dependence term), we have by Lemma I2.BI that 

[ i^'^\N,) fi^^\dx) = [ (EA)(iV,) (EA)(dx) = (EA)'(iV(A')) (3.6) 
Jx Jx 

and 

e(^Jj^{N,) Udx)^ =E(4'Ha^W)) =E(A[2l(Ar(^))). (3.7) 
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For the third integral (weak dependence term) some more work is necessary. The term that 
we would like to estimate is 



2 ; sup 



c 



a{dx), (3.8) 



where u'-'^'' is the density of /x^'^'. Equations (lA.Sh and (IA.4j) from the appendix imply that, for 
almost every x £ X and for C S M{N^), 



fNl{g) PiMds) = / E(Q^'(aUc) I e = a)f{a) Pr{da) 



and 



(3.9) 



4"uw(^? + '^-) PiMdQ) = / ^7T{x)Q^^\a\Nc) I e = a + 5,)f{a + 5,) Pi{da), (3.10) 

where Q[?\fT) = Eeca,£,ec(nsee^'(^)) (nsG<7\£,(l - every a e m* and every p e E, so 

that [(p, cr) 1-^ (5|^?'((t)] is £^ A/'-measurable. By Equation p.Sp we have furthermore that 



It 



(TT) 



(x) = / + Pi,Ns{d0) = f E(^(x) \i = a + 5^)f{a + 5^) Pi{da), (3.11) 



using that Q|^(^c)('7) = 1 for every a G Oa*(iV^). 

The absolute value term in (13.81) can then be estimated as 



[fNU^}(0 + - /S'(^?)^*"nx)] Pl,iVc(df?) 



/ E{TT{x)Q^^\a\Ns) U = <T + 5.)/(a + 4) Pi(da) 

- / E{Q^^\a\Ns)\( = ^)f(.<^) Piida) ■ I H<x) \ ^ = a + 5^) f{a + 5^) Pi{da] 



< 



+ 



E{7r{x)Q^^\a\N^) \ ^ = a + 5,) /{a + 5^) Pi{da) 

- / E(Q(7'(aUc) \ ^ = a + 5,) E{tt{x) \ i = a + 6,)f{a + 6^) Pi{da 

[ E{Q'^\a\Ns) \ ^ = <^ + ^x) E{7t{x) \ C = a + d.,)fia + 4) Pi{da) 

- / E{Q^^\a\N^)\^ = a)f{a) Piida) ■ [ E{tt{x) \ ^ = a + 6^) f{a + 6,) Pi{da 

[ coY{7T{x),Q^^\a\Ns) U = + 5.)/(a + 5^) Pi{da 

^{Qc\0cxt\Ns) I ^Uext = Qext) E(vr(x) I ^Ui„t = ffint + (^x) 

/AextUAi„t(£'ext + Qint + S^) Pl,Ai^,{dgint) A,Aext (c^£'ext) 

E(Qc^(£»extU|) I CUext = £'ext)/Aext(£'ext) A.Aextlc^e'ext) 

?n(Aoxt) 

/ E(7r(x) I e|Ai„, = 0nt + 5x)/A„,(f?int + <5x) Pl,A„,(d0nt) , (3.12) 



+ 
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where Condition ()3.2p was used for the last equahty. Note that Q^'(cr|Arc) depends on vr only 
via vtIjvc, and does so in a i-(A''^)-measurable way, where £{N^) = a(^[p i— > p{x)]; x G A'^^) is the 
canonical a-algebra on E{N^) := {p\Nf, P £ E}. 

The first summand on the right hand side of Inequality (|3.12p can then be bounded further 

as 

/ cov(7r(x), Q^T^aUi) \ ^ = a + 6,) f (a + 6^) Pi{da) 

cov{tt{x), h{Tr\N^) I ^ = fj + 4.) f{a + 5^) Pi{da) 

cov(7r(x), /i(7r|Arc) \ ^ = a + 5x) f{(J + S^) Pi{da) 



< I ess sup 

Im h:E{Ni)^[0,l] 



ess sup 

^ h:EiN-)^{0,l} 

where the essential suprema are taken over all £'(A'^^)-measurable functions with values in [0, 1] 
and {0, 1}, respectively. The third line is obtained by 

|cov(x,y)| < ||y||L^|cov(x,sign(y))| 

< 2||y|lL^maxf|cov(X,I[y > 0]) L |cov(X,I[y < 0]' 



for all random variables X G Li and Y G L^o, and for Y := E(X [ 1") — EX (see [TO], Section 1.2, 
proof of Lemma 3), where we set X := tt(x) and Y := /i(vr|Arc) — 1/2. 

For the second summand on the right hand side of Inequality ()3.12p . we use the notation 
FciQext) ■■= ^{Qc\Qext\Ns) \ CUext = S'ext) and G{gint + Sx) ■■= E(7r(x) | ^Ui^t = ftnt + Sx) , and 
bound it as 



Jcxt , 



< / G{0int + Sx) sup 



1 



G{Qint + Sx) 



f{Qcxt,Qmt +Sx)\ Pl,Ac^t{dQcxt) Pl,Aintid0mt), 



Pl,Ai,,id0i, 



where the supremum is taken over 7\A(Acxt)-iiieasurable functions. The third line is shown by 
setting Fo(£)cxt) := l[/(£'ext, ftnt + Sx) > O] and noting that Fq : ^{Ac^t) [0, 1] is measurable 
and maximizes the absolute value term in the second line. 
Thus the total estimate for the weak dependence term is 



2 / sup 



c 



a{dx) 



< 2 



7(x) a{dx) + / (3{x) a{dx). (3.13) 



Pluging ([321), ([32]), and (IXT3D into Theorem [2£] yields statement (i), and, since = EA(Af), 
also statement (ii). □ 



Proof of CoroUary \3.E[ We aim at applying Theorem 13. Al for R = 0. Clearly Condition (13. ip 
holds for any neighborhood structure. By the independence of ^ and vr we have ^(vr) as a 
regular conditional distribution of tt given the value of ^ and we see that Condition ()3.2p is 
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satisfied, that = E(7r(x))0(x) for almost every x, and that Inequality ()3.4p simplifies to 
(|3.5|) by Equation (j2.4p . Using the representation ^ = X^l^i from the proof of Lemma [2.61 
we have furthermore that 



EA{A) = E(E(^Tr{Si)l[S, G A] 



E( ^E(7r(5,) I e A] 

E7r(x) ^i{dx) 



for every A G B, and by the analogous computations that 



EA[2](B)= / E{tt{x)tt{x)) fi,2]{d{x,x)) 
Jb 

for every B G B^. Based on these results we can apply Theorem I3.AI and obtain the upper 
bounds stated. □ 



Proof of CoroHarv \3.Cl We apply Corollarv lS.BI for the point process ^ := which has hered- 
itary density / := fx with respect to Po(Leb^|;i:') and expectation measure jl = miLeb^lx, 
where all of these objects are interpreted as living on X (as opposed to living on y and be- 
ing trivial on 3^ \ ^). Consider as neighborhood structure {Nx)xex given by := fl X, 
write N!^ for the complement of Nx in X, and set N{X) := {{x,y) G Af x 3^; y G N^} and 
N{X) := G X^; y G Nx}^ which are measurable by the fact that the Nx are translated 

copies of a single measurable set. Denoting the conditional density based on / by g, we obtain 
for the 0(x)-term 



IE|ff(x;f|^,) - mi 



\fNsu{x}i0 + hciQ) mi\ P^^^cidg) 



\f{x\N:,)u{x}{Q + Sx) - fx\N^iQ)mi\ Pi^xxnMb) 



m(M^\iXUN^)) 



{f{M^\N^)U{x}{B + B + ^x) 



< 



- f{M^\N^){Q + q) mi) Pi,Ai,\(xuN^){dQ) 



Pl,X\NjdQ) 



^9'{x;C\m,\nJ - mil, 



(3.14) 
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and thus by Corollary I3.BI that 

dTv{^{^p\x),MP^iLeb^W)) 

<p^ f fii{N^) fii{dx)+p'^fi[2]{N{X))+p [ E\g{x;i\j;^,)-mi\dx + 
J X Jx 

<p^ pi{N^) fiiidx) +p'^p[2]{N{X)) +p I E\g'{x;^\M^\Nj dx. 
Jx Jx 

Statement (i) follows from this by noting that /M[2](iV(A:')) = 'ml\X\}C{N) (see Equation ()4.5p ) 
and using the various spatial homogeneities that were required. Statement (ii) is obtained 
likewise, using additionally that EA(,^) = pmi\X\. □ 

Proof of Corollary \3.M From the definition it is clear that the total variation metric is not 
affected by changes of scale of the state space, so that 

dTv{^{Cnf^T^),Po{p^^^K^^)) =dTv{-^{i.),Po{p^^')). (3.15) 

The definition of di and the inequality required for d^ imply that di{QiK7p^ , Q2K7p^) ^ di{gi, Q2) 
for all Qi, g2 & whence, by Equation (12. 3p . 

d2{^{Cn^T')^Po{p^^^K^')) <d2{^{(^),Po{p^^^)). (3.16) 

With Equations (13.150 and ()3.16|) it is seen that all the results from Theorem l3.Al to Corollary l3.CI 
remain correct if we do the proposed replacements; in particular, the upper bounds stated follow 
immediately from Theorem I3.A[ □ 



4 Applications 

We study two applications for a fairly general point process ^ here. The first one concerns the 
thinning of ^ by covering it with an independent Boolean model. This is up to a few technical 
adjustments the setting that was used in Section 3.3 of [TS]. We present this application in order 
to illustrate to what degree the results of the current article improve on the main distance bounds 
in [18], and give new insight into the high intensity limit behavior. The second application deals 
with a Matern type I thinning of ^. We present it as an example where the rather involved 
/3-term is non-zero and can be reasonably simplified. The bound is compared to a result in 
where the same thinning was considered for the special case that ^ is a Poisson process. 

In this whole section we consider a metric do on that is generated by a norm, and use 
notation of the form B(x,r) for closed do-balls in and M'^{x,r) for their complements. Write 
furthermore ]B;f(x, r) := B(x, r)n*Y for the corresponding balls in Af and B^(x, r) := X\Mx{x,r). 
We call the subset X of admissible if it is compact, of positive volume, and has a boundary 
dX that is of volume zero. 



4.1 Thinning by covering with an independent Boolean model 

The details for this situation are as follows. 

Model Setting 1. Suppose that X C is admissible and that ^ is a point process on X which 
has a density / with respect to Pi := 'Po{Leh^\x) and finite expectation measure n = ni with 
density u. Let q G [0, 1], and take H to be a stationary Boolean model (see [H], Section 3.1) 
on R-^ whose grains are do-balls of random but essentially bounded radius, denoting by /i > 
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the intensity of the germ process and by Ri € L^o the radii of the grains (which are i.i.d.). This 
means that H takes the form 



[jBiYi,Ri), (4.1) 



i=l 

where Yi are the points of a Po(/iLeb'^)-process that is independent of (i?j)jgN- Assume fur- 
thermore that ^ and H are independent, and define a retention field by 7t{lo,x) := ql[x "(i^)] 
for uj Q and x (z X. Thinning with respect to vr corresponds to deleting all the points that are 
covered by H, while retaining uncovered points independently of one another with probability q. 



We aim at applying Corollary 13. Bl in this setting. Assume without loss of generality that 
P[i?i > 0] > (otherwise Proposition I4.AI below is easily checked directly), and remove from 
H(lij) n X any lower-dimensional parts, stemming either from balls with radius zero or from 
balls that only just touch X from the outside, by taking the closure of its interior in R^. 
Note that this does not alter the distribution of the obtained thinning, because only a set of 
volume zero is removed in this way and because and vr are independent. As a consequence of 
Proposition EiHiv), where y = X and S = Q^n;f , we obtain then that vr has an evaluable path 
space. Let Nx ■= Mx{x,f) for some f > 2\\Ri\\l^ and every x € X, which implies independence 
of tt{x) and tt\n^ and hence that we can choose 7 = in Inequality ()3.5p . We set furthermore 
r := so that = K[R^). By the fact that the capacity functional T3 of the Boolean 

model H is given by 

Ts{C) :=P[HnC/ 0] = 1 -exp(-/iE(Leb^(B(0,i?i) + C))) 

for any compact set C C (see [21] . Equation (3.1.2)), we obtain for the expectations in the 
upper bound of Corollary 13. Bl 

E7t{x) = q{l - Ts{{x})) = qe-l^nmRi)\ (4.2) 

and 

E(7r(x)7r(x)) =g2(i_TH({x,£})) =(7Vii^l«(°'^i)^«(^-^'^^)l. (4.3) 



As earlier, we use absolute value bars for a measurable subset of M to denote its Lebesgue 
mass. Defining az) := \M{0,1)\ and 6 : ^ [0,1] hy b{y) := E\M{0, Ri)\M{y, Ri)\/E\M{0, Ri)\, 
we then have the following result. 

Proposition 4. A. Under Model Setting 1 laid down above and letting r > 2||i?i||i^, r := 
WRiWlb! ^'^^ Nf:{X) := {{x,x) E X'^; dQ{x,x) < f}, we obtain that 



dTv{^{^.),Po{l^n) < q'e-^'^^^'%l{N,{X)) 



e 



'li {l+b{x—x))a]jr 



D 



H[2] {d{x,x)) 



where /z''^' = qe '^"^^^/ii and 



:= sup sup 

xex ceJ\f{B^{x,f)) 



c 



[fM^{x,f)u{x}{Q + Sx) - fM^{x,f){Q)u{x)] Pifi^(x,f){dg) 



(4.4) 
□ 
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Remark 4.B. Under Assumption 1 made for Proposition 3.G in [18], the above estimate can 
be bounded by const- (f^lA'Ig^e"'^"^^^ + \X\qe~^^'^°^° (3^^^^^^ , which makes it somewhat better 
than the one in Proposition 3.G (if we accept P^^^"^^ as a natural substitute for /3(™'^^(m) in 
[TS] and apply Equation (j3.15p ). also since the result is formulated in the stronger (irv'-™etric 
instead of d2- 

However, the main point worth noting here is that the derivation above is considerably simpler 
and more elegant than the one for Proposition 3.6, because we do not have to worry about 
covering discretization cuboids. For the same reason we are easily able to work with balls that 
are based on other metrics than the Euclidean one and can write down the explicit constants in 
the upper bound. 

If we assume that ^ is second order stationary (i.e. the restriction to X oi a second order 
stationary point process C on M^), the rather complex second factorial moment measure can 
be replaced by a term involving the corresponding reduced moment measure. Second order 
stationarity means that the second moment measure /i2 of ( is locally finite {^2{B) < oo for every 
bounded measurable B C M^) and invariant under translations along the diagonal {{x,x); x G 
M^} (see [8], Definition 10.4.1), and implies stationarity of the expectation measure, so that 
^1 = miLeb^ for some mi G It follows from Lemma 10. 4. Ill in [S] that there is a measure 
JC on M.^ (unique if mi > 0) such that 



h{x,x) fj.i2](d{x,x)) = mf / / h{x , x + y) }C{dy) dx (4.5) 
for every measurable function h : M+. Hence 



< mllX] [ e-'it^+^^J')!"^"'' }C{dy). (4.6) 



i,f) 

If C is stationary and mi > 0, it can be seen by Equation (j4.10p (see [21] , beginning of Section 4.5) 
that the measure /C is given as 

JC(B) = — ECo(S) (4.7) 
mi 

for every Borel set B C M^, where Co denotes the reduced Palm process of C given a point in 
(see [13], Lemma 10.2 and Section 12.3), so that K(q{B) can be interpreted as the expected 
number of points of C in -B given there is a point in 0. The measure /C is usually referred to as 
second reduced moment measure, although some authors prefer defining it as mi (or even mf) 
times the above measure. Set furthermore K{f) := /C(B(0,f)) for every f S which, if do is 
the Euclidean metric, defines Ripley's iC-function. 

We examine the situation of Corollary 3.H in [18], waiving two technical conditions that were 
needed there, but insisting on second-order stationarity in order to bring the second summand 
in the upper bound in a nicer form. 

Model Setting 1' . Suppose that J C MP is admissible, that n G N, and that X = = 
K~^(J), where Kn{x) = {l/n)x for every x G MP . Let ^ be a second order stationary point 
process on X which has density / with respect to Pi = Po(Leb^|;f) and expectation measure 
ji = Hi = miLeh^\x for some mi G M+. We assume that ^ is the restriction to X of one and 
the same point process C, on for every n, and suppress the index n in any quantities that 
depend on n only by virtue of this restriction. Choose a sequence {qn)n& with 1/n^ < 'Zn < 1, 
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and a sequence (H„)„gN of stationary Boolean models on R-^ of do-balls with radii Rn,i € Lc 
(i.i.d. for every n) and germ process intensity Zi > such that 

1 , . n.\l/B 



= (7^1og(g„n^)) ' . (4.8) 

Assume that ^ and H„ are independent for every n G N, and define retention fields by 7r,i(a;, x) := 
qnX\x for a; € and x ^ X. Let furthermore f„ > 2||i?„^i||i^, and note that 

/^Wft-i = = miLeb^\j. 

By Inequality gj]) and as dry (^(e.^/tn^), Po(miLeb^|j)) = ^tv (^(^.J, Po(^Leb^|Af)) 
by Equation p.lSp , we have the following consequence of Proposition I4.A[ 

Corollary 4.C. Under Model Setting 1' we obtain that 
dTy(^(C.„K;'),Po(miLeb^|j)) 

< ml\J\aDr-Y + m\\JW^ [ (g„n^)-''(^) /C(dy) + 2| Jj/^t"^ 



<\J\{mlaD{ — ) +mlqnK{rn) + 2p, 



n 



D 



(sup) \ 



where ^'f^^^ was defined in Equation (14. 4p . □ 

Remark 4.D. Note that Assumption lb) in [18] implies that K{fn) = 0{f^) for n — > oo. Com- 
pared with the corresponding result in [18j (Corollary 3.H) we therefore have again somewhat 
better bounds with explicit constants that were obtained in a more direct way. 

A rather nice result can be derived from Corollarv I4.CI in the Poisson case. 

Corollary 4.E. Under Model Setting 1' and the additional assumptions that C, is a Poisson 
process with expectation measure miLeb^ , that qn = 1 for every n € N, that do is the Euclidean 
metric, and that ||-Rn,i||Loc = 0{rn) for n — > oo, we have 

dTy(^(e.„'^^'),Po(miLeb^!j)) = 0((log n)-(^-i)) for n ^ oo. 

Remark 4.F. The following (partly heuristical) arguments suggest that the order claimed in 
Corollarv 14. El is sharp for mi > 0. Assume for simplicity that J = [0, 1]^. 

For = 1 it is readily understood that ^7r„K~^ cannot converge in distribution to a Poisson 
process. The reason is that the uncovered part of M in the domain of the contraction, i.e. M\l^i^, 
is made up of intervals whose lengths are exponentially distributed with mean l/li (no matter 
how the Rn^i are distributed), so that the probability of having two or more points within the 
first uncovered interval that lies completely in does not depend on n. Hence we have a 
constant positive probability that the first two points of the thinned contracted point process in 
are within distance l/(Zin), say, which cannot be true for a sequence that converges towards 
a homogeneous Poisson process. 

For D >2 the situation is more complicated. By Theorem 1 in [11] (compare also statement (ii) 
on page 244) , it can be seen that the uncovered "chinks" of H„ in the domain of the contraction 
have a volume that is of order 1 / (log(n^)) ^ for large n, so that the argument of the constant- 
sized chinks is not valid for general D. Heuristically, the order of the chink volumes together with 
the fact that the number of chinks in a bounded measurable set is Poisson distributed (see |llj . 
p. 244, statement (i); note the slightly different scaling) suggest that we can think of the process 
^■K„K-n^ for n large as a compound Poisson process YlJ=i '^j'"^'^^?") '^ith intensity of the Poisson 
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process X^i^i '^gC") °f order (log(n-^))^ ^ and i.i.d. clump sizes Z^^\ for which P[z|"^ > 1] is 

of order (log(n^)) and P[zj"^ > 2] is of order (log(n^)) ^'•'^ (by the fact that ^ is a 

Poisson process). It is easily seen that such a process converges towards a homogeneous Poisson 
process ?? as n ^ oo by noting that c^tv ^^s("))>-^(E^=i ^ > ll-^^w)) ^ 0, 
and YliJ^i — Vi but that its convergence rate in the total variation metric is 

bounded from below by ¥[3i G {1, . . . , Vn} '■ > 2], which is of order (log(n^)) or, 
what is the same, order (logn) 



-(D-l) 



Proof of Corollary \4-I^ Our starting point is the first upper bound in Corollary 14. CI where we 
set fn '■= 2||i?n,i IIloo • Since C is Poisson, the third summand is zero, whereas the first summand 
is clearly 0((log 71.)"^^"^^) . We investigate the integral in the second summand. We have 
IC = Leb^ by Equation ()4.7p in combination with ^(Cd) = -^(C) (see [H], Proposition C.2). 
Define b : [0,2] [0,1] by b{u) = |]B(0, 1) \ ]B(y, 1)|/|B(0, 1)|, where y is an arbitrary element 
of with \y\ = u, and note that b{y) > 6(2|y|/f„) for y £ M{0, Tjj). Since do is the Euclidean 
metric, it can be shown that there is a constant k > such that b{u) > for every u G [0, 2]. 
Writing uj d for the surface area of the unit sphere in , we then can bound the required integral 
as 

/ (gnn^)-^(^) lC{dy) < [ (n^)-^^y\/rr. Leb^(dy) 

JM{0,r„) 7B(0,f„) 

- l'\n^rr^-' dr 



K 







< -D (log(n )) (1 



where the last inequality follows from multiple integration by parts. Since r„ = 0((log(n^))^/^) , 
we thus obtain that also the second summand in the first upper bound of Corollary 14. CI is 
0((logn)-(^-i)). □ 



4.2 Thinning by Matern type I competition 

Again we base our retention field on a random closed set S, but this time we choose a situation 
where H is completely determined by the point process ^. The resulting procedure is the one 
used for the construction of the Matern type I hard core process, in which a point is deleted 
whenever there is any other point within a fixed distance r. The details are as follows. 

Model Setting 2. Suppose that r > and that X,X' C MP are two compact sets, where 
X is admissible and M{X,r) C X' . Let furthermore ^ be a stationary point process on X' 
(i.e. the restriction of a stationary point process C, on MP) which has density / with respect to 
Pi := Po(Leb^|;i") and a finite expectation measure fi = fii = miLeh^\x' for some mi G M+. 
By defining ^ on X', but considering the thinned point process only on X, we avoid boundary 
effects, which would lead to more complicated notation because of spatial inhomogeneities in 
the thinned process. 

In order to have a H whose realizations are closed sets that are jointly separable, we proceed 
as follows. Write ^ as YlJ=i^Si, where V and Si are cr(^)-measurable random elements with 
values in and X', respectively, and denote by Tj the do-distance between Si and its nearest 
neighbor. Let then 

Bi{Lo) := {y G X'; ^ min(rj(w), r) < do{y,Si{io)) < r} 
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for every u; G $7, and set 

oo 
i=l 

Choose q E [0, 1] and define a retention field on X' by setting 7r(u;,x) := ql[x if u; € O 

and X £ X and 7r(ix', x) := if w E and x £ X' \ X . Note that 

7r(cj,s) =gl[s E Af,^(a;)(]B(s,r)) = O] (4.9) 

for CO £ ^} and s E '^(w), where B(x,r) := B(x,r) \ {x} for every x E Hence, on X, thinning 
with respect to vr corresponds to deleting all those points that see any other point of the process 
within distance r (regardless whether this point is itself deleted or not), while retaining points 
that do not have this property independently of one another with probability q. 

This time, we aim at applying Theorem 13. Al for the state space X' . By Proposition I A.El f iv) . 
vr has an evaluable path space (after removing from 'B{uj)r)X possible lower-dimensional parts by 
taking the closure of its interior in , which has no influence on the distribution of the resulting 
thinning). Since vr is completely determined by ^ we have the corresponding Dirac measure as 
a regular conditional distribution of vr given the value of ^. Condition ()3.2p is satisfied for a 
catchment radius R = r, so that := Mx'{x,f) for some f > 2r is a legitimate choice for 
the neighborhoods of strong dependence. We can furthermore choose 7 = in Inequality (13. 4p . 

Write for the reduced Palm process of ^ given a point in x, and ^ for the second-order 
reduced Palm process of given points in x and x (see [13], Section 12.3, pp. 109 &: 110; note 
that = H[n] for obtaining the distributions of the n-th order reduced Palm processes). The 
first and second order Campbell-Mecke equations state that 

h{x,i-5x)i{dx)^ = j m{x,il) iiiidx) (4-10) 

and 

e( f h{x,x,i-5x-5x)i^^\d{x,x))]= [ E/i(x,x,e-,i)/U[2](d(x,x)) (4.11) 

for non-negative measurable functions h. These equations follow immediately by standard ex- 
tension arguments from the defining equations of Palm processes (see e.g. [T^, Equation (C.4), 
for the first one). We then obtain by Equation ()4.9I) that 

EA(yl) = e( /" gl[^(B(x,r)) = O] C{dx) 
xJAnx / 

= Q f F[cl{M{x,r)) =0] fii{dx)=miq{l-G{r))\AnX\ (4.12) 
JAnx 

for any A£B = B{X'), where G : M+ ^ [0, 1], G{f) := P[C(i(B(0,f)) > l] = P[Ci(B(x,f)) > l] 
for arbitrary x E M"^, denotes the nearest neighbor function of C,, i.e. the distribution function 
of the distance from a "typical point" to its nearest neighbor, which is a frequently used tool in 
spatial statistics; see e.g. ^ (Section 3.4), or [Hj. In a very similar way, using in addition 
Equation (|4.5p to obtain the last equality, we have with Nf{X') := {(x,x) E X'^; dQ{x,x) < f} 
and Nf{X) := Nf{X') Ci X'^ that 

EA^^^NfiX')) = e( [ qH[C{M{x,r)) = 0]l[e(B(i, r)) = O] C^^\dix,i)) 

= q^ I l[(io(x,x) >r]P[e-,s(B(x,r)UB(x,r)) = O] /i[2] (d(x, x)) 
JNf(X) 



m 



jq^ / / (1 - G2,y{r)) }C{dy) dx, (4.13) 



X J {M{0,f)\B{0,r))n{X-x) 
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where G2,y : M+ ^ [0,1], G2,y{f) := P[Cd,j,(B(0, f) U B(y, f)) > l] = P[Ci,,+,,(B(x, f ) U ]B(x + 
2/)^)) ^ l] for arbitrary x € M^, are the natural two-point analogs of the G-function, with 
y G M^. 

Finally, the term P{x) = f3f{x) is zero for x ^ X' \ X , and can be evaluated for x ^ X as 



-.{x) = [ 



gl[ftnt = 0] 



?t(B(x,r)) 

,-D 



[/(fext, ftnt + 4)| -Pi 



-Our' 



^,(x,f-r)) 



^,{x,r~r)) 

fTB{x,r)\M'y,{x,f-r){^x',£,\lB'i^,{x,f-r)) - fM{x,r){^x) 



(4.14) 



where we set /Ai„t | A.^t (^'mt; &xt) := /AextUAi„t (fet + £'int)//Aext (^'ext) if /Aext(£'cxt) > and 

/Aint I Aext (ftnt; e'ext) := otherwise. 

The expectation in the last line of ()4.14p is much simpler than the general expression we 
have for /3, but is typically still hard to estimate. A more directly applicable estimate, which, 
however, is very rough and works only with point processes that are not too extreme in a certain 
sense, is given as follows. Choose f := 2r, and assume that / > Pi-almost surely. We then 
obtain with r/cxt ~ Pi,M'',{x,r) that 



(3{x) 



q^Qint = 0] 



\f{Qext, Qint + Sx)\ Plfip{x,r){dQext) Pl,-B{x,r){dQi 



int) 



01 



x,r)) J'yi(M^,(x,r)) 
Ik|B(x,r) = 0] f{cr\M^,(x,r) + ^x) - /b^, (x.r) (o"Ib^, {a;,r))/B(a;,r) (<Ja;) Pl{da) 
f^" ,{x,t){'^W ,{x,r))fv,{x,r){^x) 



= q [ I[cT(B(x,r)) =0] 
< mig(l — G'(r)) 



1 



f{(^\B^,{x,r) + ^x) 
fm'i, {x,r){'>]ext) fM{x,r){^x) 



/(??ext + Sx) 

for x € X, where we used in the last line that 



(4.15) 



miE(/i(x,e-)) = / h{x,a)f{a + 6x) Piida) 



for almost every x (z X' and every non-negative measurable function h, which is a consequence 
of the first Campbell-Mecke equation and of Equation (j2.7p with = {x}. We assume that 
the ||-||i;,^-term is bounded by a constant M G R+ (that depends neither on r nor on x). Two 
examples where this is satisfied are given at the end of this subsection. 

Pluging ()4.12p to (j4.15p into Theorem l3.A( i) and choosing now f := 2r everywhere yields 
the following result. 

Proposition 4.G. Under Model Setting 2 laid down above, we obtain that 
dTy(if(C.),Po(/x(-')) <m2|;i^|2^aflrV(l-G(r))' + m?|A'|(72/' (i _ ^2,^,(0) IC{dy) 

JB{0,2r)\B(0,r) 



Ix 



fM{x,r)\MS^,{x,r){^x'-,^\M^,{x,r)) ~ fB{x,r){^x) 



dx, 
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where /x*'^' = mig(l — G{r))Leb \x- If in addition / > Pi-almost surely and the IMIloo'^^'""^ 
in Inequality (j4.15p is uniformly bounded by M G M+, then the last summand can be estimated 
further by 

mi\X\q{l-G{r))M. ^ 

Remark 4.H. In order to obtain an integrand that does not depend on x in the last summand 
of the above bound, we can either proceed as in Corohary 13. CI applying Inequality ()3.14p and 
as a consequence replace B^,(x,r) by Mx \IB(x,r) for bounded "outer neighborhoods" 
that are shifted copies of one another and that all contain the set X'] or we can proceed as 
in Remark 13. Dl using a Gibbs construction on the whole of MP and as a consequence replace 
B^,(x,r) by \]B(x,r). 

If ^ is a homogeneous Poisson process and q = \, then is the usual Matern type I hard 
core process restricted to X , and the above bound takes especially simple form. 

Corollary 4.1. Under Model Setting 2 and the additional assumptions that Q is a Poisson 
process and q = 1, we have 

dTy(^(e.),Po(/iLeb^U)) < WaDr''ll + nil\X\ [ g-m,lB(o,OuBfo,r)| 

iB(0,2T-)\B(0,r) 

where li := mie"''"^"^''^ . 

Proof. The first inequality follows directly from Propositon 14. Gl by the fact that ^ is a Poisson 
process, and hence the last summand is zero and the one- and two-point G-functions can be 
easily computed by -S^(Cd) = ^{C) and ^(Cd,y) = ^{C) (see [U\, Proposition C.2, for the 
one-point case; the two-point case is proved in the analogous way). 

The second inequality is a consequence of |B(0,r) UB(y,r)| > ||B(0,r)| for |y| > r. The 
latter is due to the fact that all the do-balls of fixed radius are translated copies of one another 
which are convex and symmetric with respect to their centers, and can be seen as follows. The 
symmetry implies that any hyperplane through the origin divides the volume of B(0,r) in half, 
while the convexity implies the existence of a supporting hyperplane Hx at every point x of the 
boundary of B(0,r), which means that Hx contains x and that B(0, r) lies completely in one of 
the closed half-spaces defined by Hx- Thus Hx — x and Hx divide into three parts, each of 
which contains half of B(0,r) or B(x,r), whence we obtain that |B(0,r) UB(a;,r)| > ||B(0,r)| 
for \x\=r. Clearly, |B(0, r) U B(y, r) | > |B(0, r) U M{x, r) | if |y | > |x| . □ 

Remark 4. J. In Theorem 6.6 of ^2] a situation very similar to the one in Corollary 14.11 for the 

special case that we choose do to be the Euclidean metric was considered. The only substantial 
difference is that in ^22j the Poisson process ^ is defined on X instead of the superset X' , which 
leads to less competition near the boundary of X and consequently to a non-stationary thinned 
process. However, this difference enters the upper bounds in [22] only insofar as balls are always 
restricted to X instead of being balls in . 

Disregarding these boundary effects, we see that the estimates in Corollarv 14 .11 are slightly better 
than the one for the total variation in [22], because our second estimate above is bounded by 
2\X\miaD{'^r)^li, which is exactly the estimate in if we adapt it to our notation. 
The main reason for formulating Corollary 14.11 however, was not to improve on this earlier 
bound, but to demonstrate that Proposition 14. G|. which holds for a much greater class of point 
processes, provides a reasonable estimate in the special case of a Poisson process. 
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We end this subsection by giving two examples of point processes for which the additional 
boundedness condition in Proposition 14. Gl is satisfied. 

Example 1. Consider a point process density of the form f{g) = X^^^f^g) for a function / that 
is bounded and bounded away from zero. One particular instance of such a density is given as 



f{g):=nX'^'(l min do{s,s)<r 



+ 7I 



min do{s, s) > f 



for every G 91 = where 7, f, A > are parameters and k > is a normalizing constant. 

Note that for 7 = 1 we obtain the density of the Po(ALeb^|;f')-process. 

Example 2 (Strauss process). Consider the Strauss process with range f G [0,r] and interaction 
parameter 7 G (0, 1] (see [Qj, Section 6.2.2). This process has a density given by 

fig) := KX^^h''^^^ 

for every g & where Cf{g) := J2s s€g syts'^ido{s, s) < f] counts the pairs of points that lie 
within distance f of one another, A > is an intensity parameter, and k > is a normalizing 
constant. Then, for a; G and cJext ^ x' '''")) ^ 



/(o-cxt + ^x) 
1 

^^l^cxt 1+1 ^Cf (Ccxt) 



< 



<n{M(x,r)) 



^ \ kcxtl + lftntUcf ((Text) p , Jrin- ^ 



^^l&xtl+l p_^. 



{x,r) 



(dgext) 



j.^,{x,r)) 

^kintl ^Itotl P^^^c r^^^^{dg^^t) Pl,M{x,r)id0int) 

x,r)) J%-B^,{x,r)) 

= K [ AI-^I Pi{da) = Ke^~\ 
J<yt 

where we used for the inequality that Cf{a + g) > Cf{(T) for all £», a G and Cf{a + 6x) = Cf{a) 
for (T G O'l and x G X with a{Mx'{x,f)) = 0. Thus, M may be chosen to be max(l, Ke^~^ — 1). 



Appendix 

In what follows we formulate and prove some of the more technical results needed in the main 
part of this article. 

A.l Density formulae used for Theorem 12. CI 

Proposition A. A. For a point process ^ on X with density f with respect to Pi and finite 
expectation measure, and for a neighborhood structure {Nx)x£X that satisfies Condition ()2.6p . 
we have 

El/ h{x,(,\Ns) Cidx)] = / h{x,0)fNsu{x}{Q + Sx) Pi,Ns{dg) a{dx) 
\Jx J Jx J'niN^) 

for every non-negative or hounded measurable function /i : ^ x — > M. 
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Corollary A.B (Generalized Nguyen-Zessin formula on compact spaces). Suppose that the 
conditions of Proposition \A.A\ hold and that f is hereditary. We then have 



for every non-negative or bounded measurable function h : X x ^ ^ M, where g is given in 
Equation (|Z!9]) . 



Remark A.C. Note that the statement of Corollary lA.BI is wrong in the case = {x} for 
all X G ^t" if / is not hereditary and ^ 0] > 0. As a counterexample consider the process 
that scatters a fixed number n > 1 of points uniformly over X (cf. Example 14.2(a)). This 
process has a density given by f{Q) = e°^'^^ ~ hence satisfies g{x, g) = for all 

X £ X and g € which makes the right hand side in Corollary lA.BI zero for every function h, 
whereas, with h{x, g) = 1, the left hand side is equal to E^(Af) > 0. 

Since the corollary does not hold generally, its use in the proofs of Theorem 2.4 and 3.6 of [2j 
and in the proof of Theorem 2.3 of ^ is not justified unless an additional condition (such as 
hereditarity) is imposed. 

Proof of Proposition \A.Al We proof the statement for non-negative h; the statement for bounded 
h follows in the usual way by decomposing h into its positive and negative parts. The Slivnyak- 
Mecke theorem (see [13], Theorem 3.2, for X C M^, which clearly can be generalized to arbitrary 
compact metric spaces) states that, for ~ Pi and measurable /i : A" x — > M+, 

E^y" h{x,r] - 6^) r]{dx)^ = j m{x,r]) a{dx). (A.l) 
Hence, setting h{x,a) := h{x.,a\Ng)f{(y + ^x), we obtain 



h{x,g)fN-,\j{x}{Q + ^x) Pi,N^{dg) a{dx) = I I h{x,a\N^)f{a + 5x)Pi{da)a{dx) 

h{x,a\Nf)f{a) a{dx) Pi{da) 



=¥.{^j^h{x,i\Ng)i{dx)y 

Proof of Corollary \A.E\ The statement follows immediately from Proposition lA.Al using that 
/Af=u{x}(^ + ^x) = g{x, g)fN^{g) for every x ^ X and every g £ OT(iV^). □ 

A. 2 Density of the thinned process 

Let the point process ^ and the random field tt be as for the definition of the thinning in 
Subsection 12.21 We assume additionally, as in Section [3l that all the realizations of vr lie in an 
evaluable path space E C [0, 1]^ and that there is a regular conditional distribution of vr given 
the value of ^ (see Appendix IA.3P . It is essential for the construction below that we use the 
same such distribution throughout (i.e. without changing it in between on P^^^-null sets), but 
insignificant, of course, which one we use. 
Set then 



q{g\a):=E(j[7T{s) J] (l - ^s)) ? = ^) 
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for almost every o" E and for g C cr. It can be easily seen that the mapping [9^^ x E ^ 
[0, 1], (q, g,p) I— > Y\g^gP{s) nse^(''^ ~P{^))] is A/"^ (8) ^'-measurable, whence we obtain that q{g \ a) 
is well-defined and that if : —>■ [0, 1], {g, g) i— > q{g \ g + g) is measurable. 

Lemma A.D. A density of the thinned process with respect to Pi is given by 

f(-){g) :=e"W / qig\g + g)fig + g)Pi{dg) 

for almost every g £ 

Proof. The well-definedness and the measurability of f^^^ follow from the measurability of ip 
defined above. 

Consider two independent Poisson processes r/,f/ ~ Pi. Take furthermore x ~ ^2 and 
Ist Xi/2 be a thinning of x with retention function p = 1/2, which corresponds to picking a 
subset of the points of x uniformly at random. Note that {rj, fj) = (xi/2) X \ X1/2) (see e.g. [H], 
Proposition 3.7, for X C M^; the proof can easily be adapted for general compact metric spaces). 

Integration of the proposed density over an arbitrary set C € A^, using Lemma 12. Al for the 
fifth line, yields 

/ f'^\g)Piidg) = e"^'^^ [ [ q{g\ g + ~g)f{g + ~g) Pi{d~g) Pi{dg) 
= e"(^)E(^I[7/ G Cjqiri \ rj + ^)/(7? + fj)" 
= e"WE(l[xi/2GC]<?(xi/2lx)/(x)' 



^"^^^/ :i^T.^[Q^C]q{g\a)f{a)P2{da) 

Y,nQ(^C]q{g\a)f{a) Pi{da) 

E{Q^^\a)\C = a)fia)Pi{da) 
= [ ¥[C^eC\( = a]fia)Piida) 

= P[e. GC7], (A.2) 

where Q^^tr) := J2eca,eec{Us€gP(^)) {Usea\oi'^ - Pi^))) every a em* and every p£E, so 
that [{p, a) ^ Qc^cr)] is £:(g)A/'-measurable and Qc\0 = £ C | ^, vr]. From Equation fO]) 
the claim follows. □ 

The proof above yields that 

/ f''\Q)Pi{dQ)= [ E{Q'^\a)\^ = a)f{a)Pi{da) 
Jc Jm 
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for every C E AA, and hence, by Equation ()2.5p . that more generally, with A G B, 

, fI\Q) PiAdQ) =11 f^^\Q+~Q) PiMds) A,A(df?) 
Ic JcJm{A<:) 



= [ r\a)P^{da) 
Jc 

^Q'~\<T)\i = cj)f{a)P^{da) 



for every C S A/'(^), where C := {(T € DT; € C}. For the last equality we used 
Y.ec.\AAlis&,Pis)){\[s&(.Uc)\s{'^ - P{s))) = 1. 



^Q'c\'y\A) I i = (T)f{<j) Plida) (A.3) 

that 

The various computations in ()A.2p remain correct (after the obvious minor modifications) if 
we add an extra point to g, yielding in the unrestricted case 

[ f^^KQ + 5x)Pi{dQ)= I y^l[Q(^C\q{Q + 5^\a + 5^)f{a + 5^) Pi{da) 
Jc J^g^^ 

= [ E{TT{x)Q^^\a)\^ = a + 6,)f{a + 5,)Pi{da) 

for every C G AA, which holds for a-almost every x ^ X. Hence we obtain in a very similar 
fashion as in Equation ()A.3p that 

J^fAu{.}(e + S.) PiAidQ) = I^E{7r{x)Q^^\a\A) \( = a + 5,)f{a + 5,) Pi{d<j) (A.4) 

for every A ^ B and every C G J\f{A), which holds for a-almost every x ^ X. 

A.3 Technical conditions on ^ and tt: evaluable path space and regular con- 
ditional distribution 

Consider a locally compact, second countable Hausdorff space y that is equipped with its Borel 
(7-algebra B = B{y). This is the most common type of space on which general point processes are 
defined. Any such space is separable, and a complete metric d can be introduced that generates 
its topology. With regard to the main part of this article, y is usually just our compact state 
space X, but it is sometimes useful to consider a natural superset of X (e.g. \i X (Z MP). 
For sets of functions y [0, 1], we introduce the concept of (locally) evaluable path spaces. 

Definition. Let E C [0, 1]-^ and let £ be the canonical cr-algebra on which is generated by 
the evaluation mappings ^'^^ : -E — > [0, 1], p i-^ p{x)-, where x ^ y. For f7 G ;S set furthermore 
E{U) := {p\u'-, P G E} and write £{U) for the corresponding cr-algebra generated by '^u,x ■ 
E{U) ^ [0, 1], p I— > j5(x), where x G U. 

(i) We call E an evaluable path space if the mapping <I> : E x y ^ [0,1], (p,x) i-^ pix) is 
£ ® ;B-measurable. 

(ii) We call E a locally evaluable path space if the mapping ^jj : E{U)xU — > [0,1], {p, x) i— > p{x) 
is £{U) ® i3[/-measurable for every U C y that is open and relatively compact. 

It can be easily seen that every locally evaluable path space is evaluable. 
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For the main results of this article we assume that vr takes values in an evaluable path space 
and that there exists a regular conditional distribution of vr given the value of ^. Neither of 
these assumptions presents a serious restriction, because they are both naturally satisfied in 
many practical applications, and if they are not, we can modify vr accordingly (provided it is 
measurable in the sense of Subsection 12. 2p without changing the distribution of the resulting 
thinning. To see this let Rhe a third of the minimal interpoint distance in which is positive 
except on a null set, and let n{io, x) := it{u), S{u))) if there is a point S{io) of ^(w) within distance 
R{uj) of X and tt{uj,x) := otherwise. We have as path space E for vr the space of all functions 
p : X ^ [0,1] that are zero except on finitely many non-overlapping closed balls of positive 
radius, on each of which they are constant. By Proposition lA.El fiii) below it can be seen, using 
the separability of y, that this is an evaluable path space. A regular conditional distribution of 
vf given the value of ^ can then be defined in a very straightforward manner, using the regular 
conditional distribution of (7r(s))sgo- given ^ = (t. 

Since the above construction looks rather artificial in many situations, we provide a few 
manageable conditions under each of which a path space is (locally) evaluable, and hereby 
substantiate the statement that an evaluable path space is naturally present in many practical 
applications. The proposition below is essentially the "path space version" of Proposition A.D 
in [18]. Where it was conveniently possible, we have generalized the conditions from to the 
space y. 

Definition. We call a set E C [0, 1]-^ separable from above [or below] if there exists a countable 
set Ti C y such that for every p £ E, every open d-ball B <Z y and every y £ M we have that 
p{x) > y for all x € n S implies p{x) > y for all x € i? [or p{x) < y for all x G i? fl S implies 
p{x) < y for all x £ B, respectively]. 

Proposition A.E. A set E C [0, l]'^ is a locally evaluable path space if it satisfies any one of 
the following conditions. 

(i) y = and there is a closed convex cone A C of positive volume such that every 
p £ E is continuous from A (see the definition in U8^ . Appendix A. 3); 

(a) Every p £ E is lower semicontinuous, and E is separable from above; 

(Hi) Every p £ E is upper semicontinuous, and E is separable from below; 

(iv) Every p £ E is the indicator of a closed subset ofy, and the family C of these closed subsets 
is separable in the sense that the subsets are jointly separable, i.e. there is a countable set 
ay such that C = C H S for every C £ C, where the bar denotes topological closure. 

Proof, (i) This follows with some minor adaptations from Proposition A.D in [18] by letting 
ft := E and defining vr to be the identity mapping on 17, so that tt{uj, x) = uj{x) for every uj £ Q 
and every x £ R^. Proposition A.D(i) in \18l yields then that : E{R) x ^ [0, 1], (p, x) i-^ 
p{x) is £{R) ® Sij- measurable for every bounded open rectangle R C MP (it is easy to see that 
the result from [18] carries over from to M^). The same result for a general bounded open 
set U instead of R follows by writing [/ as a countable union of open rectangles. 

(ii) We can essentially reproduce the prove of Proposition A.D(ii) in [1^. Let y £ [0, 1) and 
U d y he open and relatively compact. We show that $^^((y,oo)) £ £{U) ®Bu. Choose a 
separant S as in the definition of the separability from above and set Q := {]B°(x, 1/n) C U;x £ 
S',n £ N}, where B°(x,l/n) denotes the open d-ball with center x and radius 1/n, and S' is 
an arbitrary countable dense subset of y. Noting that p~-^((y,oo)) n C/ is open by the lower 
semicontinuity p £ E, we have 




29 



(iii) Apply (ii) to E' := {1 - p ; p € E}. 

(iv) We apply (iii). It is evident that every p € E is upper semicontinuous. Separability of E 
from below is inferred from the separability of C as follows. First note that the definition has to 
be checked only for y = 1, because p ^ E takes only the values and 1. Let E be a separant for 
C as in statement (iv), and let B be an open ball in y. Then, for C G C and p = Iq-, Pi^) < 1 
for all re G n S implies that 5 PI S fl C = 0, hence B C {C Ci S)'^. Since B is open, this implies 
B C int((Cn S)'') = [CnT^y = C^, where mt{A) denotes the interior of the set A for any 
Acy. Thus p{x) < 1 for ah x e B. □ 
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